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ABSTRACT 

This is the first in a series of papers that introduces a new paradigm for understanding the jet in M87: a 
collimated relativistic flow in which strong magnetic fields play a dominant dynamical role. Here we focus 
on the flow downstream of HST-l — an essentially stationary flaring feature that ejects trails of superlumi- 
nal components. We propose that these components are quad relativistic magnetohydrodynamic shock fronts 
(forward/reverse fast and slow modes) in a narrow jet with a helically twisted magnetic structure. And we 
demonstrate the properties of such shocks with simple one-dimensional numerical simulations. Quasi-periodic 
ejections of similar component trails may be responsible for the M87 jet substructures observed further down- 
stream on 10 2-3 pc scales. This new paradigm requires the assimilation of some new concepts into the astro- 
physical jet community, particularly the behavior of slow/fast-mode waves/shocks and of current-driven helical 
kink instabilities. However, the prospects of these ideas applying to a large number of other jet systems may 
make this worth the effort. 

Subject headings: galaxies:individual: M87 — galaxies: active — galaxies: jets — methods: numerical — 
MHD 



1. INTRODUCTION 

M87 is a nearby giant elliptical galaxy (Virgo A, 3C 274, 
NGC 4486 ) located at the cent er of the X-ray-luminous Virgo 
cluster (Fabric ant et al. 111 980). host of the first extragalactic 
jet discovered dCurtis I ll 9 1 8h . The one-sided jet emerging 
from the nucleus of M87, where a 3.2 (±0.9) x 10 9 solar mass 
black holfl resides (Mac chetto et al. I fl997). has been well- 
studied on a wide range of wavelengths from radio to X-rays 
dOwen et aL]ll989tlBiretta et al I [1991 IPerlman et al. | | 1999h 
Biretta et all 1999): lJunor et al. lll999HMarshall et al. fbooj 
Wilson & Yang|2002l:lPerlman & Wilson 120051: lH arris et al. I 
20061: iLv et al. 1 120071 : iKovalev et al. 1 120071) . Because of its 
proximity (D = 16 Mpc: lTonry||199lT) . which gives a lin- 
ear scale of 78 pc arcsec -1 , the M87 jet is one of the best 
candidates to investigate relativistic outflows in extragalactic 
systems. 

One of the most remarkable features of the M87 jet is 
the innermost bright knot G, lying about 1" f rom the core 
in the Very Large Array (VLA) observations dOwen et al. I 
[l989h . That region has been resolved by Hubble Space 
Te lescope (HST) into a structured complex known as HST- 
1 dBiretta et al. 1 1 19991) . It is located around 0.8-1" (pro- 
jected) from the core (or ~ 260 — 320 pc de-projected for 
a viewing angle of ~ 14°; Wang & Zhou 2009). It ap- 
pears almost stationary, but component velocities downstream 
of the complex are highly relativistic, with a range 4c — 6c 
dBiretta etal. 1119991) . Indeed, the component of HST-1 that is 
furthest upstream (i.e., HST-ld) is stationary to within the er- 
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rors (< 0.25 c), and has been identified as the poi nt of origin 
of the superluminal ejections dCheung et al. 120071) . As of this 
writing, no other observations have detected any superluminal 
components upstream of the HST-1 complex. Therefore, the 
observations obtained up to now, paint a surprising picture: no 
evidence for highly relativistic velocities between the core and 
HST-1, a stationary knot at HST-ld, and then, suddenly, su- 
perluminal motion of components immediately downstream 
of HST-1. 

The structure of the jet downstream of HST-1 (1 - 18" or 
0.1 - 1.5 kpc in projected distance) can be characterized by 
trailing clumps or knots of bright gas {HST-1 to C) with an ap- 
parent deceleration t o sub luminal speeds (around 6c to 0.5c) 
dBiretta et al. JT_995. 1999f) and filamentary structures ("wig- 
gles/kinks") Jowenetal.l [l989t ISparks et al. I [[996). The 



high overpressure in the s ynchrotron gas and the highly po- 
larize d helical filaments dOwen et al. 1 119891 IPerlman et al. 
1999) indicate the existence of underlying ordered magnetic 
fields with a 3-dimensional helix seen in projection; magnetic 
fields, therefore, appear to play a role in determining the M87 
jet structure even on large scales. The magnetic field vec- 
tors in the knots HST-1, D, A, and C are perpendicular to the 
jet direction, indicating the presence of longitudinal compres- 
sio n by a shock front and/or a tightly wou nd magnetic he- 
lix dOwen et al. II 1989b IPerlman et al. ll2003l) . Detailed broad- 
band (from radio through optical to X-ray) spectral shape of 
the knots (HST-1 to C) in the M87 jet favors the scenario in 
which synchrotron emission dominates the radiation and in 
situ particle acceleration (by the first-order Fermi process) al- 
most certainly occurs in the large scale M87 jet (b oth within 
knots and outside them) (Perlman & Wilson 1120051) . 

The purpose of the paper is to investigate the dynamics 
downstream of HST-1 using relativistic magnetohydrodynam- 
ics (MHD). We introduce here the concept that the superlu- 
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minal components in M87, and perhaps extragalactic jets in 
general, are the result of MHD shocks produced in helically 
twisted, magnetized relativistic outflows. Our exclusive focus 
here on the flow just downstream of HST-l, a unique dynam- 
ical behavior of a pair of sub/super-luminal knots, constitutes 
the first in a series of papers on the MHD paradigm for the 
entire M87 jet. The paper is organized as follows. We outline 
our model in fj2] In fj3]our numerical method is introduced, 
and in S|4] we show our numerical results. Discussions and 
conclusions are given in |J5] 

2. MODEL DESCRIPTION 

We propose an MHD model of the M87 jet that extends be- 
yond, or downstream, of HST- 1 that resolved by the VLA and 
HST observations (1 — 18" or 0.1 — 1.5 kpc in projection). 
The structure of the jet downstream of HST-l can be char- 
acterized by a conical shape z on r (where r is a jet radius 
and z is a dista nce from the nucleus) wit h an o pening angle 
of 0q„en ~ 6° (Ide Vaucouleurs & Nieto ||1979t IOwen et al. I 
1980; Reid et al. lll982tlOwen et al. II1989I) . The high degrees 
of radio and optical polarization in b oth the knot (40% — 6 0%) 
and interknot (20% — 40%) regions dPerlman et al. Ifl99l in- 
dicate the presence of coherent magnetic fields (highly or- 
dered) on large scales associated with the underlying jet. 

Under the assumption of the minimum energy condition, 
the kn ots themselves ap pear to be significantly overpres- 
sured dOwen et al. 1 1 19891) w ith respect to the ambient ther- 
mal gas ([Yo ung et al71 l2002h . but the interknot regions do not 
(Spar ks et al. 1119961) . In order to maintain a conical stream- 
line of the adiabatic jet within a uniform ambient gas (an iso- 
thermal King profile with a core radius R c ~ 18"), the fields 
may have to be much stronger and more highly ordered (i.e., 
a force-free configuration with a 3-dimensional helix) than a 
weak and tangled field at the equipartition level (~ a few of 
100 /iG). Magnetic fields, therefore, appear to play a crucial 
role in determining the structure of the M87 jet even on larger 
scales over 100 pc. 

In the present paper, we do not examine the entire three 
dimensional structure of the M87 jet, described in Sjl] this 
will be treated in forthcoming papers. Instead, we focus on 
the unique properties of the ejected superluminal components 
(ICheunget al. 1120071) by using a simple one dimensional ap- 
proach. ICheung et al. I (12007) reported that some time be- 
tween 2005 December and 2006 February, one of the newly 
ejected components, called HST-lc, split into two bright fea- 
tures: a faster moving component (cl: 4.3c ± 0.7c) and a 
slower moving (c2: 0.47c ± 0.39c). If the bright knot is a 
shock, then such a pair of separated components naturally can 
be identified as forward/reverse modes. We here model this 
interesting feature by relativistic MHD simulations. In addi- 
tion, we suggest that a similar expansion might be occurring 
in earlier HST observatio ns (HST-le: 6.00c ± 0.48c / HST-l 
East: 0.84c ± 0.11c) by lBiretta et al. I (119991) . 

Furthermore, the Knot D complex located just downstream 
of the HST-l complex, also displays distinct features (DE - 
DM - DW) with different propagation s peeds and a sidew ays 
oscillation of the brightness ridge-line (Birett aet al. 11995b . It 
has been pointed out that four quasi-periodic knot complexes 
downstream of H ST- 1 (called D, E, F, and I) are at roughly 
uniform spacing (Biretta & Meisenheimer 1993). Thus, re- 
current ejections of superluminal components at the HST- Id 
site may happen with a certain quasi-periodic time scale. We 
suggest that these bright components can be considered as 
similar paired structures to 7/5T-lc2/cl and HST-l East/e: DE 
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FIG. 1. — Schematic view of our model of the 87 jet in the VLA scale. 
The structure of the jet downstream of HST-l can be characterized by trailing 
MHD quad shocks that may be intermittently generated at the HST-\ com- 
plex. The system of MHD quad shocks consists of the forward fast (FF) and 
slow (FS) followed by the reverse slow (RS) and fast (RF) shocks. 

-DM - DW, E - EF F - 1, a nd A - B - C further dow nstream 
[refer to lOwen et aU (119891) and lBiretta et all d 19991) for def- 
initions of the know labels]. 

The stationary knot in HST-l d may be identi- 
fied as the recollimation shock dStawarz et al. I 120061 ; 
Bromberg & Levinson 2009). Once the recollimation shock 
is formed at a finite distance, the reflection of supersonic 
flows may emit a shock component downstream of the 
recollimation point. A Steady model of relativistic MHD 
mu lti-layer outflows (rel ativistic jet and non-relativistic wind) 
by iGracia et al. I (|2009) nicely reproduces su ch observed 
prope rties of M87 as an asymptotic collimation dJunor et al. I 
1999) with a feature brightened at around 100 pc (projected 
distance) by an over-collimating MHD flow. 

Of particular intere st is the brigh t est em issi on structure A - 
B -C. As reported bv lOwen et al. I (11989ft and lPerlman et al. I 
( 1999), knots A and C have certain similarities: i) bright trans- 
verse linear features (normal to the jet axis) indicative of a 
shock front (Biretta e t al. 1119831) : ii) dominance of transverse 
magnetic field suggesting ordered helical magnetic compo- 
nents. Visible side-to-side oscillation is also observed be- 
tween these knots, and magnetic vectors appear to follow the 
fluctuating jet axis in this part (including knot B). 

We point out a sudden enhancement of emission at knot A, 
at the upstream edge of that knot, indicating a reverse shock 
and a rapid drop in emission at the downstream edge of knot 
C, suggesting that it is the correspon ding forward shock (see, 
Figure 8a, 10a of lOwen et al. I [19891) . Particle acceleration 
is associated with both the forward and reverse modes. In 
addition, the substructure between the upstream edge of knot 
A and the downstream edge of knot C indicates there may be 
another pair of (forward/reverse) shocks or rarefaction waves. 
The downstream edge of knot A or knot B may be a reverse 
mode feature, while the upstream of knot C may be a forward 
mode feature. 

Sudden changes in magnetic vector orientation strongly im- 
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ply the existence of MHD fast/slow mode waves. The trans- 
verse component of the magnetic field B (B±) increases 
across a fast mode front (the normal component B\\ re- 
mains unchanged through an MHD oblique shock and thus 
the magnetic pressure p m increases), while B± decreases 
across a slow mode front (p m decreases). Thus, we inter- 
pret the brightest emission structure A - B - C to be a trail 
of quad relativistic MHD shocks (forward fast/slow and re- 
verse slow/fast) (perhaps, two slow modes can be rarefaction 
waves) generated in a helically twisted super-fast magnetic 
flow. Such a flow can be generally expected at further down- 
stream from the AGN in the MHD jet theory and emissions 
of MHD shocks from the stationary knot feature has been 
inspected by MH D simulations of strongly magnetized jets 
dLind etal. 111989b . The presence of four shocks, instead of 
two, is due to the presence of B^, which serves a role equiva- 
lent to B± in the simpler planar oblique shock case. 

The transverse field component is dominant in some knots 
(HST-1, D, A, and C), indicating strong longitudinal compres- 
sions by a passing shock. In the past, a transverse orienta- 
tion of the magnetic field was thought to be due to longitudi- 
nal compression of roughly random (tangled) magneti c fields 
(weak ly polarized) by a hydrodynamic-like shock (Laing 
1980) (i.e., this is es sentially a perpendicu lar shock). How- 
ever, as is shown in iPerlman et al. I (1 1999b . the high degree 
of polarization is confirmed in interknot regions of the M87 
jet on scales of 10 2-3 pc. So the classical picture of a weak, 
random jet magnetic field may be in conflict with these ob- 
servations. It is natural to conclude, therefore, that underly- 
ing magnetic fields in both the interknot and knot regions are 
systematically ordered with helical (longitudinal + azimuthal) 
components (three-dimensional helical structure). Longitudi- 
nal compression of a helical magnetic field also will produce 
an enhanced azimuthal (i.e., projected transverse) component. 

High energy particle acceleration by the first-order Fermi 
process is generally believed to occur in extragalactic jets 
(Blan dford & Ostriker 1 119781) . A relativistic particle energy 
distribution [n(E) oc E~ s ] steeper than S — 2 would be 
needed to produce the radio - optical - X-ray synchrotron 
spectrum in the M87 jet. Synchrotron model fits from radio 
through optical to X-ray data produce <5 = 2.2 at all ene rgies 
and all locations along the jet ( Per lman & Wi lson 20051) and 
about S = 2.36 on average dLiu&Shen 112007b . These agree 
very well with the conditions n eeded for a diffusive sh ock ac- 
celeration (DSA) [5 = 2- 2.5: lKirk & Dendy I(l2001b l. 

The geometry of a perpendicular MHD shock (this is a 
so-called "magneto-acoustic shock", a particular case of the 
fast-mode MHD oblique shock with only magnetic compo- 
nent along the shock surface B±) is not suitable for a DSA; 
negligible p article acceleration has been confirmed b y PIC 
simulations (lLangdon et al. Ill 9881: iGallant et al.l fl992). Fur- 
thermore, in the standard picture of Fermi acceleration, the 
DSA does not work for slow-mode MHD shocks because the 
magnetic field strength decreases across them, and therefore 
they cannot act as magnetic mirrors for the upstream parti- 
cles. Thus, only fast-mode MHD oblique shocks are probably 
responsible for the observed non-thermal emissions. 

In summary, we suggest that each pair of bright knots on 
10 2-3 pc scales in the M87 jet is a quasi-periodic event pro- 
duced in the stationary HST-id knot as a consequence of rec- 
ollimation processes of a converging super-fast magnetosonic 
jet. They are trailing quad MHD shock wave fronts (forward 
fast/slow modes and reverse slow/fast modes), which propa- 



gate as super/sub-luminal components in a highly magnetized 
and twisted, relativistic ou tflow powered by a non-linear tor- 
sional Alfven wave train dMeier et al. Il200l1) . A schematic 
view of our M87 model is shown in FigureU 

The physics of MHD shocks in helically twisted magne- 
tized flows are still not very well known in the astrophysi- 
cal jet community. Despite the fact that jet dynamics is in- 
herently three-dimensional, much of the physics of these jets 
can be obtained from simple 1-D simulations of flow along 
a cylindrical shell. Therefore, in this paper we in vestigate 
observed properties of HST-lc (Cheung et al. II2007T) using a 
simple one-dimensional approach that suffices in illustrating 
the basic principles. We demonstrate that observed proper 
motions of forward (superluminal) and reverse (subluminal) 
knots can be reproduced precisely in our relativistic MHD 
simulation model. 

3. NUMERICAL METHOD 



Based on observ ations (see, e.g., lOwen et aT~l 1 19891: 
IPerlman etaTT l999). we assume that the magnetic field plays 
a fundamental role in determining the flow properties of the 
M87 jet all the way up to scales beyond 100 pc in projection. 
We impose axisymmetry and choose a cylindrical coordinate 
system (r, </>, z) whose axis coincides with the symmetry- 
axis. We model the dynamical behavior of observed HST- 
lc knot which split into two distin ct features as sub/super- 
luminal knots (ICheung et al. 1 12007). We here consider initial 
phase of the knot ejection and separation (propagation of in- 
dividual knots with constant speeds) and thus dynamics of 
the flow can be described in so-called 1.5 -dimensional ap- 
proximation along a cylindrical shell of radius ro (which is 
assumed to be rigid), which allows the quantities to vary in 
the z-direction, and which also allows for the influence of az- 
imuthal effects. 

Our system obey s ideal, special relativistic MHD 
dLichnerowicz 1 1 1 967 ) that consists of the baryon mass and 
energy-momentum conservation laws in the absence of a 
gravitational field, Maxwell equations in CGS-Gaussian units, 
and Ohm's law: 



dQ dF 

dt dz 



0. 



IP 



^hV^/c 2 - E r B z /{ATTc) 
Q= 7 2 hV z /c 2 + E r B^/(4Trc) 
B 

fh-p+(E 2 r +B 2 +B 2 M{^) 

JPVz 

^hVtpVz/c 2 - B^B z /(4tt) 
F(Q) ee ^hV^/c 2 +p+ (E 2 + B 2 - B 2 z )/(8n) 

cE r 

j 2 hV z +cE r B^/(4ir) 



(D 



(2) 



(3) 



The above equations denote the mass, the azimuthal (<fi) mo- 
mentum, the axial (z) momentum, the azimuthal induction, 
and the energy equations respectively. 

Here p and p are the rest-mass density and gas pressure in 
the fluid frame. V = (V^, V z ) is the fluid velocity and 7 
is the associated Lorentz factor 7 = 1/(1 — V 2 / c 2 ) 1 / 2 . c 
is the speed of light, and T = 5/3 is the specific heat ratio. 
B = (£?0, B z ) denotes magnetic field as measured in the 
laboratory (galaxy) frame. (Note that B z is constant by V • 
B = and never changes in time.) h = pc 2 + Tp/(T — 1) 
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is the relativistic enthalpy, and E r [= (V^,B Z — V z B$)/c] is 
the radial component of the electric field as measured in the 
laboratory frame. 

We normalize physical quantities by unit length scales Lq, 
density p , velocity Vq in the system, and other quantities de- 
rived from their combinations, e.g., time as Lq/Vq, etc. This 
normalization does not change the form of the basic equa- 
tions. A factor of An has been absorbed into the scaling 
for both E r and B. The resulting set of time-dependent, 
fully conservative equations for Special Relativistic MHD 
(SRMHD) is solved by a finite volume method (FVM). We 
use a newly designed, hybrid flux (HF) scheme (Nakamura 
2010) that consists of a unique hybridization of a Godunov- 
type and centered difference-based fluxes to achieve low-level 
numerical diffusion. (Typical 1-D Riemann problems exhibit 

^HF/^HLL _ q g an( j jHF^Roc _ q 7> where g j g ^ relative 

error dete rmined from an L\ norm). MUSCL-TVD recon- 
struction (I van Leerlll979|) and MUSCL -Hancock predictor- 
corrector time marching ( van Leer 1984) schemes, with a van 
Albada limiter (Ivan A lbada et al. 1982), are implemented to 
maintain second-order accuracy in both space and time. 

4. NUMERICAL RESULTS 

We adopt ro = 1 pc as the unit length Lq (and as the radius 
of HST-ld) and c = 3x 10 10 cm s _1 as unit velocity Vq. An 
electron density n e — 0.17 cm~ 3 and a temperature kT = 
0.8 keV at the Bondi radius - 120 pc (lAllen et al. Il2006h 
are taken as central ISM properties. These correspond to a 
density pisu{= fJ-"m p n e ) = 1.7 x 10~ 25 g cm~ 3 (m p is the 
proton mass and p = 0.61 is the molecular weight for full 
ionization) and temperature Tism = 9.3 x 10 6 K. We also 
assume a unit density po = 1.7 x 10~ 27 g cm -3 , which is 
two orders smaller than pisu (light jets are generally believed 
to be present in radio galaxies and quasars) (see, e.g., iKrause I 
120031) . This gives p c 2 = 1.5 X 10~ 6 dyn cm -2 as the unit 
pressure. The unit time to becomes ro/c = 10 8 s = 3.2 yr, 
and the unit magnetic field Bq is (Airpoc 2 ) 1 / 2 = 4.3 mG. 

The computational domain z e [—0.04, 2.0] (pc in a di- 
mensional scale) is resolved with 5100 grid points. A Rie- 
mann problem, which possesses two uniform initial states (1: 
left and r: right) separated at z = 0.0, is considered to investi- 
gate the evolution of the system. Our fiducial model consists 
of a super-fast magnetosonic, relativistic jet inflow (7 ~ 6.22) 
into a nearly force-free, weakly twisted {B^jB z ~ 0.33) 
medium flowing with a sub-relativistic speed (7 ~ 1.07). It 
produces a set of quad relativistic MHD shocks. The follow- 
ing initial conditions are prescribed: 

(p, V*, V z , B+, B z , pf = (1.0, 0, 0.987, 0.8, 2.4, 0.256) 

on the left-hand side (-0.04 < z < 0.0) and 

(p, V4,, V z , B4, B z , p) 1 = (1.0, 0, 0.36, 0.8, 2.4, 0.256) 

on the right-hand side. Time integration, using a CFL number 
of 0.8, is followed until t = 2.0 (~ 6.4 yr) to examine the 
earliest phases of relativistic MHD shock propagations. 

The above choice of initial values has been inferred care- 
fully from various observational co nstraints. By assumin g a 
viewing angle 6» v ~ 14° at HST-l dWang & Zhoull2009l) . a 
maximum pattern speed of a faster moving component HST- 
7c 1 (including an error) c an be estimated from its apparent 
speed /3 app = 4.3c ± 0.7c dCheung et al. 112007b as 

pattern = ~ R ^ . fl * 0.982, (4) 

Papp cos o v + sin V 



Contact discon. (CD) 




FIG. 2. — Space-time (z, t) diagram of logarithm of the proper density ■yp 
(vertical axis), as measured in the laboratory (galaxy) frame. Quad MHD 
shocks; FF, FS, RS, and RF, and a contact discontinuity (CD) are labeled. 



where j3 — V/c. It would seem far more natural to have the 
pattern speed tied to the jet fluid speed, with /3 pa ttorn % Pfluid 
dBirettaetaL~lll995l) . We note pattern = 0.989 corresponds 
to 6 V = 19°, an upper limit of the viewing angle for pos- 
sible solutions to have the sup erluminal motion /3 a DD = 6.0 
as seen in HST observations (Bir etta et al. 1119991) . On the 
other hand, an ambient motion in the vicinity of HST- 1 com- 
plex can be constrai ned by the stationa ry feature of HST-ld 
with /3 app < 0.25 dCheung et al. 112007b as ppattcm < 0.516. 
The magnetic field vectors in projection around HST-l are 
perpendicular to the jet direction and the timescale of opti- 
cal and X-ray variab ility of HST-l requires \B\ ~ 10 mG 
dPerlman et al. Il2003h . Taking into account a viewing angle, 
the constraint for the proje cted .B-vectors to b e perpendicular 
to the jet can be expressed dAsada et al. 1120081) as 

> sin 6 V ~ 0.24. (5) 

B z 

Thus, our simulation model is reasonably guided by observa- 
tions. In the following, we examine our results. 

Figure|2]shows the propagation of MHD wave fronts in the 
proper density jp. Quad MHD shocks and a contact discon- 
tinuity (CD or entropy wave), all with constant speeds, are 
clearly visible. Relative to a reference frame that co-moves, 
and co-rotates, with the jet plasma near the CD, these waves 
propagate in both the forward (F) and reverse (R) directions. 
Here we adopt the convention of counting shocks beginning 
with the one farthest from the origin of the disturbances (HST- 
1). Two of the four shocks, the first and the fourth, are forward 
fast-mode (FF) and reverse fast-mode (RF) shocks, respec- 
tively. The other two, the second and the third, are forward 
slow-mode (FS) and reverse slow-mode (RS) shocks. 

Snapshots of various quantities at t = 2.0 are illustrated in 
Figure [3] As seen in (a) and (b), the gas is compressed across 
the first (FF) and second (FS) shocks. In crossing the third 
shock (RS) the gas is expanded by a smaller ratio than the 
FS, while it is expanded much more strongly in crossing the 
last (RF) shock rather than the FF as clearly shown in (b) [it 
may not be visible in (a), but this is due to the frame of the 
measurement; in the rest frame of fluid elements, the distri- 
bution of p, instead of jp, has a similar tendency with p]. As 
a result, the gas accumulates in the region between the sec- 
ond and third shocks. As one moves from large to small z, 7 
increases with gradual steps throughout the first, second, and 
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third shocks, but largely increases at the last shock towards 
the injection level 7 ~ 6.22 shown in (c). Similarly, from 
(d), increases across the first shock, decreases across the 
second one, increases again across the third shock, and finally 
decreases across the fourth one. From (e), changes as well, 
in a way consistent with the increased twist between the first 
and the second shocks, and reversed between the third and 
fourth shocks Q 

We define the plasma-/? (a ratio of the gas pressure to the 
magnetic pressure) in the rest frame of the fluid element: 



/3 P 



2p 



B 2 



Bih 2 



(6) 



This enables us to compare the importance of the magnetic 
forces with the plasma forces in the proper rest frame of the 
fluid. As shown in (f), the unshocked region further down- 
stream of the FF is highly magnetically dominated /3 p ~ 0.08, 
which is prescribed as the initial condition. f3 p increases 
slightly in crossing the FF (both the gas and magnetic pres- 
sures are enhanced similarly). However the accumulated gas 
region (RS-CD-FS) is heated strongly (the gas pressure in- 
creases, while the magnetic pressure decreases) and thus at- 
tains near equipartition (f3 p ~ 1.4). At the forward part of 
the shock quad (FS-FF) (3 p is still very small (~ 0.1), but at 
the reverse part (RF-RS) it has increased by nearly one or- 
der (~ 0.9). In the region behind the shock quad, after the 
four MHD shocks have passed, f3 p decreases to a low level 
that is prescribed in the initial condition. Given the fact that 
the relevant physical quantity is the fluid frame value of the 
plasma-/3 (/3 P : eq. |6), and that B^/ B z < 1, f3 p is not appre- 
ciably different from that in the initial values adopted in the 
galaxy frame. Note that the relativistic length-contraction ef- 
fect produces a weak compression of the structures between 
FS-FF rather than between RF-RS when viewed in the labo- 
ratory frame (which is also the computational grid on which 
the simulation is done). 

Strengths and propagation speeds of the four shocks re- 
main constant with distance as they propagate in our coor- 
dinate system: axial propagation (z-direction) in a uniform 
medium (constant sound and Alfven speeds) in a fixed-radius 
cylindrical shell. Individual speeds are estimated as Vff ~ 
0.97 c, Vfs ~ 0.94 c, Vrs ~ 0.73 c, and V" RF ~ 0.67 c, 
re spectively. For a vi ewing angle of 9 ~ 14° at HST- 
1 dWang& Zhou 1 120091) . the faster component HST-lcl has 
~ 0.97 c, while the slower component HST-lcl has ~ 0.67 c. 
As is mentioned in a separation of observed super/sub- 
luminal components can be identified as distinct proper mo- 
tions of two fast-mode MHD shocks, Vff and Vrf- There- 
fore, our n umerical model is r easonably consistent with ob- 
servations (ICheung et al. Il2007h at a quantitative level. 

We performed a simulation with an underlying velocity (the 
left state in the Riemann problem we solve) that was con- 
strained by the maximum velocity of the HST-ld (< 0.25c) 
allowed by the observations. There is no evidence for an ad- 
ditional underlying flow other than the velocity of the HST-ld 
component itself. That is, the flow that produces the superlu- 
minal motion of HST-lcl is the flow, and that is our model: 
the ejection of a new MHD jet from HST-1 at a speed much 
greater than any underlying flow is what prod uces the recent 
observations of HST-1 complex suggested by Cheung et al. 

4 Note that the region FF - FS and the region RS - RF are counter-rotating 
when viewed from a frame that rotates with the plasma near the CD. 




FIG. 3. — (a)-(f) log(7p), log(p), 7, Vj>, B^, and /3 P (plasma-,3), respec- 
tively, shown at t = 2.0. Only the region 1.2 < z < 2.0 is displayed. 
Note that panels (a)-(e) are measured in the rest frame of the galaxy, but the 
panel (f) is measured in the rest frame of fluid element. Each discontinuity is 
labeled in (a) (see also Fig. [JJ. 



(120071) . The FF shock that mainly determines the Lorentz fac- 
tor for the HST-1 complex is HST-lcl. What we learn from 
the observations is that the RF shock {HST-lcl) is sublumi- 
nal; a property predicted by the quad shock model. This con- 
stitutes a constraint despite the fact that the exact speed of the 
RF shock is determined by the model parameters which have 
been specified only within a given range. However, we em- 
phasize that a very wide range of those parameters predicts 
sub-luminal motion for the RF shock. 

In a forthcoming paper, we will study the propagation of 
quad relativistic MHD shocks in a conical geometry (z oc r) 
with an increasing cross section that is compatible with the 
M87 jet on 10 2-3 pc scales. This will allow us to study decel- 
eration of the shock s as th e jet propagates on 10 2-3 pc scales 
dBiretta et al. Ifl995l[l999h . 

5. DISCUSSION AND CONCLUSIONS 

Figure [4] shows the variations in the magnetic helical lines 
of force over ~ 1 radian of the cf>-z plane as the jet (propa- 
gating in the z direction) passes through each of quad MHD 
shocks. Magnetic field strengths \B\ in the regions upstream 
of RF (©) and downstream of FF ((T)) are about 10.9 mG 
(B^/Bz ~ 0.33). As mentioned above, the azimuthal field 
components in the inter-shocked regions are amplified by 
compression, and they dominate over the axial component in 
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FIG. 4. — Vai'iation of the helical magnetic field in the (f)-z plane (B^, B z ) 
and in the frame of the galaxy, as the jet flow (in z) crosses the quad MHD 
shocks (labeled at right). (Compare the right-hand panel with Fig. Only 
~ 1 radian of the pattern is shown, which, of course, is periodic over the 
full 2ir radian circumference of the cylinder. The pitch angles of the field 
upstream of the RF and downstream of the FF (both are "shock-upstream") 
are equal (in this frame). Changes in the pitch angle are caused by jump 
conditions at the different shocks. Each state between quad MHD shocks is 
identified by numbers 



the rest frame of the galaxy; B$jB z and \B\ are ~ 2.3 and 
~ 25.7 mG (FS - FF: ©), ~ 1.2 and ~ 16.0 mG (RF - RS: 
0), respectively. On the other hand, in the rest frame of the 
fluid element, only the forward part of inter-shocked region 
(FS-FF: ©) has large magnetic pitch (B 4> /j)/B z ~ 1.3. 

One of the important insights that we derive from our nu- 
merical model is an asymmetric structure of inter-shocked 
regions such as gas compression and magnetic pitch angle 
pitch [= tan -1 {B^/j I B z )j (in the rest frame of the fluid el- 
ement) in pre-shock regions of the two fast shocks: 6> pitch ~ 
17° downstream of the FF (Q), and ~ 3° «/?stream of the RF 
((D) (#pitch in the frame of the galaxy is same at the down- 
stream of the FF and the upstream of the RF, as is shown 
in Fig. 4). Efficiencies and mechanisms of relativistic parti- 
cle acceleration depends crucially on both the magnetization 
and the magnetic obliquity of the upstream plasma. Parti- 
cle acceleration is mostly mediated by DS A for quasi-parallel 
shocks (6* pitch < 10°), but shock drift acceleration (SDA) is 
the main acceleration mechani sm for larger magnetic obliq- 
uity (ISironi & Spitkovkyll2009t) . In the rest frame of the fluid 
element, the density increases by r cmp . ~ 3.42 at the RF and 
r cmp ~ 1.90 at the FF (where 7- cmp . = P2I P\ is the shock 
compression ratio and p\, 2 are the density ahead and behind 
the shock). The slope 6 in the particle energy distribution 
for P SA is determined as S = (r cmp . + 2)/(r cmp . — 1) dBell I 
1978) and thus <5 ~ 2.23 is expected at the RF, which agrees to 
the theory and observation described in fj2] In a forthcoming 
paper, we will examine these aspects (DSA/SDA) in details to 
identify observed properties of the M87 jet in 10 2-3 pc scales 
as a context of quad relativistic MHD shock system. 

As mentioned in [21 the flow downstream of HST-l has a 
conical structure with a constant opening angle until it reaches 
knot A. Poloidal magnetic field varies as B z oc r 
region, while toroidal magnetic field varies a s B^ o c r 
With measured radii tq ~ 1 pc near HST- ld (jCheun 
2007) and n 
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33^2 - 44.1 pc at knot A dOwen et aO l98 
Sparks et al. 1 119961) . combined with the B-field strengths 
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(i?0 ~ 12.39 mG and B z ~ 10.32 mG) at the reverse feature 
(RF - RS: 0) in the simulation, we estimate B^i ~ 280.9 - 
373.2 pG and B zi ~ 5.3 — 9.4 pG as a counterpart of knot 
A. This strongly toroidall y dominated field agrees with the 
polar ization observations (lOwen et al. 1119891 iPerlman et al. I 
1999). Using the upper limits on inverse Compton radiation 
imposed by the HESS and HEGRA Cerenkov telescope ob- 
servations of the kpc scale jet t o estimate the magnet ic field 
strength in the brightest knot A, Stawar zet al. I (12005) obtain 
a "safe" lower limit of |£?| > 300 pG, indicating a departure 
from the equipartition value (from the synchrotron spectrum 
of the knot A); the magnetic field energy density within the 
brightest knot is very likely higher than the energy density 
of the radiating ultrarelativistic electrons. This implies that 
interknot regions are likely to be extremely magnetized. On 
the other hand, a reasonable upper limit to the field strength of 
knot A can be constrained by the total power Lj of the M87 jet 
from the relation 7rr 2 (V x B) x £?/(47r) < Lj/2, where Lj — 
few xlO 44 erg s" 1 dOwen et all 120001) . and we assume an 
equipartition between the matter energy and Poynting fluxes. 
In a situation in which V z ~ c > V$ V r , B^ 3> B z > B r , 
the dominant term in the axial (z) direction, the main carrier 
of electromagnetic energy, is B^V z /(Att) ~ B^c/ (Air). This 
gives l-BK-B^) < 1 mG. Thus, these constraints are consistent 
with our derivation \B\ ~ 289 — 370 pG at around knot A by 
considering an expansion of jet cross section. 

Sideways oscillations of the brightness ridge-lines down- 
stream of the HST- 1 complex are closely followe d by the fila- 
menta ry structure of distributed magnetic fields (lOwen et al. I 
1989).' Some of these filaments can be interpreted as a 3- 
dimensional helix seen in projectio n, produced by growin g 
current- driven instabilities (CDIs) (Nakam ura et al. I 2001). 
Indeed, INakamura & Meier I (120041) found that MHD shocks 
play an important role in triggering helical kink CDIs: while 
rotation of the plasma may stabilize the helical kink instabil- 
ity locally, MHD shocks (particularly in the region A - B - 
C) can rapidly alter this stabilizing rotation (Figure [3]), sud- 
denly causing the magnetized plasma to violate the Rruskal- 
Shafranov stability criterion. Further downstream of knot C, 
the knot-like features disappear, and the jet becomes diffusive 
with strong side-to-side oscillations and bending. This may 
also be a good example of the growing "external mode of the 
CDIs beyond the X-ray cluster core where a separation of the 
paths between the jet forward and return current occurs due 
to the rapidly decreasin g density of the external thermal gas" 
(INakamura et al. 120071) . We shall examine the above topics in 
forthcoming papers by treating our quad MHD shock model 
in 1-D conical and, eventually fully 3-D, flow environments. 

The idea that magnetic fields are instrumental in the forma- 
tion and propagation of jets in active galactic nuclei dates back 
four decades. Despite a recent growing consensus on this no- 
tion stemming from the results of numerical simulations of 
magnetohydrodynamic flows near black holes, the precise dy- 
namical role of magnetic fields in observed parsec and kilo- 
parsec jets has remained uncertain. In this first in a series 
of papers, we have combined SRMHD numerical simulation 
with observation to explain superluminal knot ejections from 
HST-l in the M87 jet. Indeed this is the first examination of 
superluminal motions in the extragalactic jet as a consequence 
of trailing MHD shocks in a relativistic flow that possesses 
helically twisted magnetic structures. We begin to develop 
a detailed picture for the jet in M87 that is grounded in the 
dynamics of relativistic MHD, ultimately suggesting that the 
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entire jet (from subparsec to kiloparsec scales) is an MHD 
phenomenon. Eventually our model may also be readily ap- 
plicable to other jet systems, and it may constitute the foun- 
dation of an MHD paradigm for astrophysical jets in general. 

Stimulating discussions with Colin A. Norman, Keiichi 
Asada, and Jose Gracia are gratefully acknowledged. M.N. is 
supported by the Allan C. Davis fellowship jointly awarded by 
the Department of Physics and Astronomy at Johns Hopkins 
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